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The effect of dipolar orientation with respect to the soliton plane on the physics of two-dimensional bright 
solitons in dipolar Bose-Einstein condensates is discussed. Previous studies on such a soliton involved dipoles 
either perpendicular or parallel to the condensate-plane. The tilting angle constitutes an additional tuning pa¬ 
rameter, which help us to control the in-plane anisotropy of the soliton as well as provides access to previously 
disregarded regimes of interaction parameters for soliton stability. In addition, it can be used to drive the conden¬ 
sate into phonon instability without changing its interaction parameters or trap geometry. The phonon-instability 
in a homogeneous 2D condensate of tilted dipoles always features a transient stripe pattern, which eventually 
breaks into a metastable soliton gas. Finally, we demonstrate how a dipolar BEC in a shallow trap can eventually 
be turned into a self-trapped matter wave by an adiabatic approach, involving the tuning of tilting angle. 


I. INTRODUCTION 

Last two decades have witnessed intensive investigations, 
both theoretical and experimental, on multi-dimensional soli¬ 
tons, especially in systems possessing nonlocal nonlinearity 
(NLNL) [1-3] such as Bose-Einstein condensates (BECs) of 
particles with permanent [4, 5] or induced dipole moments 
[6] (dipolar BECs), photo refractive materials [7], nematic liq¬ 
uid crystals [8, 9] and others [10-12]. Two-dimensional (2D) 
optical bright solitons [13] and three-dimensional (3D) light 
bullets [14] have been reported in non-local media. In addi¬ 
tion, nonlocal interactions significantly influence inter-soliton 
collisions [4, 7, 15, 16] and may lead to the formation of soli¬ 
ton complexes [17, 18]. 

Bright solitons in BECs have been realized experimentally 
only in quasi-lD alkali-atom gases [19-21], in which short- 
range contact-like interactions result in a cubic nonlinear¬ 
ity equivalent to self-focusing nonlinearity in a Kerr media. 
Whereas in these media higher-dimensional solitons are un¬ 
stable, 2D solitons may become stable in dipolar BECs for a 
sufficiently large dipole-dipole interaction (DDI) [4, 5]. Thus, 
the realization of dipolar BECs of chromium (Cr) [22, 23], 
dysprosium (Dy) [24] and erbium (Er) [25] opens fascinating 
perspectives for the realization of 2D BEC solitons. The re¬ 
alization of 2D solitons remains however an open challenge. 
On one hand the observation of the 2D solitons (for a dipo¬ 
lar orientation perpendicular to the soliton plane) proposed 
in Ref [4] demands the inversion of the sign of the DDI by 
means of rotating fields [26]. On the other hand, anisotropic 
solitons, proposed in Ref. [5] for a dipolar orientation on the 
soliton plane, demand no external driving, but additional chal¬ 
lenges may arise from their high anisotropy [27], in particular 
their vulnerability to collapse instability even in a strictly 2D 
scenario [28]. 

In this paper we study the effect of dipolar orientation with 
respect to the soliton plane (see Eig. I) on the properties and 
stability of 2D bright solitons in dipolar BECs. Introducing 
the tilting angle a with respect to the normal vector of the 
quasi-2D trap plane [29-31] allows for accessing previously 



Figure 1. (a) The schematic setup of dipolar BEC confined in the xy 
plane with a strong harmonic confinement along the z axis, (b) The 
dipoles are polarized in the xz plane with tilting angle a w.r.t. z axis 
and 8 is the angle between the dipole vector and the radial vector r. 


disregarded regimes of interaction parameters for soliton sta¬ 
bility as well as to manipulate the soliton anisotropy in a more 
controlled manner. These features may enhance the possibili¬ 
ties of realizing 2D BEC solitons experimentally in the state- 
of-the-art dipolar BECs. 

Solitons may be created by driving a BEC into phonon- 
instability (PI). This has been done in non-polar ID BECs 
by tuning the short-range interactions from repulsive to at¬ 
tractive [21]. In contrast, in 2D and 3D PI leads to collapse 
in non-dipolar BECs [32, 33]. On the contrary, the post PI 
dynamics in 2D dipolar BECs is characterized by the forma¬ 
tion of a transient gas of bright solitons [28]. The resulting 
solitons exhibit an intriguing dynamics that crucially depends 
on the nature of DDI on the soliton plane. Isotropic solitons 
attract each other in any direction, undergoing fusion during 
collisions and eventually becoming a single large soliton. The 
scenario is different for anisotropic solitons, as the interaction 
between the solitons is anisotropic in the plane, which makes 
them fuse only when they are colliding along the dipolar axis. 
Here we extend these studies to the tilted case. As we show, 
the post-PI dynamics is always characterized by a transient 
stripe pattern that eventually breaks into bright solitons. 

The paper is structured as follows. In Sec. II we discuss 
the model and the corresponding non-local Gross-Pitaevskii 
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equations (NLGPEs). In Sec. Ill, we examine the Bogoli- 
ubov excitations of a 2D homogeneous BEC, identifying the 
regimes of PI as a function of the interaction parameters and 
the tilting angle. Sec. IV is devoted to analyzing the stability 
and the properties of the soliton. The stable/unstable regions 
of Cr, Er and Dy BECs are discussed as a function of system 
parameters. In Sec. V we demonstrate how to prepare a 2D 
soliton by varying the tilting angle from a soliton unstable to 
a soliton stable region in the case of a Cr BEC. Einally we 
conclude in Sec. VI. 


II. MODEL 


We consider a BEC of N particles with magnetic or electric 
dipole moment d, oriented in the xz plane forming an angle a 
with the z axis, using a sufficiently large external field (Eig. 1). 
The DDI potential is Vdix) - ^^(1-3 cos^ 8)lr^, where 0 is the 
angle formed by the dipole vector d s t/(sin a x + cos a z) and 
the radial vector r, with gd ^ Nd^ being the strength of the 
dipole-potential. At low-enough temperatures the system is 
described by a non-local Gross-Pitaevskii equation (NLGPE): 
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where J (ir|'P(r, f)p = 1 and g - AnfP'aNIm is the coupling 
constant that characterizes the short-range contact interaction, 
with a the i-wave scattering length. Eurther, we assume a 
strong harmonic confinement along the z direction with a fre¬ 
quency Wj and no trapping in the xy plane, hence Vt{r) - 
mat^z^ 12. This trapping is sufficiently strong such that the sys¬ 
tem remains in the ground state, <po(z) - exp(-z^//j)/ 
of the harmonic oscillator along the z-axis with 4 = 
factorizing the BEC wave function as 'P(r) = i^(x,y)(f>Q(z). In 
that case, the physics of dipolar BEC becomes quasi-2D, with 
the restriction that \fi 2 D\ ^ where jj. 2 D is the chemical 
potential of the 2D gas. Employing this factorization, convo¬ 
lution theorem, the Eourier transform of the DDI potential. 
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and integrating over dz, we get an effective 2D NLGPE; 
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Figure 2. The PI regions (shaded) of a 2D dipolar BEC as a function 
of p and the tilting angle a (in radians) for (a) g > 0 and (b) g < 0. In 
the former case post PI dynamics is characterized by the formation 
of 2D solitons, while in the latter case the BEC becomes unstable 
against local collapses. The dashed line in Fig. 2(a) indicates the 
magic angle a„, across which the characteristics of DDI should be 
reversed in order to stabilize the 2D bright solitons. 


where we use dimensionless polar coordinates (k = 
4 yjkl + ky and 0k), and erfc(x) is the complimentary error 
function. Note that for a = 0 and a = 7t/2 Eq. (3) reduces to 
the cases discussed in Refs. [4] and [28], respectively for the 
realization of isotropic and anisotropic solitons. 


III. PHONON INSTABILITY 


In this section, we calculate the low energy spectra (Bo- 
goliubov excitations) of a 2D homogeneous dipolar BEC with 
tilted dipoles. The homogeneous solution of Eq. (3) is 
i/'(x,y, f) = ■\/n 2 D exp[-//r2Z)f/4], with 1120 the homogeneous 
2D density. The 2D chemical potential is obtained as 



where p — gdl\g\ determines the ratio between the DDI 
strength and that of the contact interactions. The elemen¬ 
tary excitations of the homogeneous solution are of the form 
dipkix, y, t) - - 1 - with the dispersion 


fk = 4wk 
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where Ek - fP' (k^ + k^ /2m. The phonon modes, e(k —> 0) oc 
yfjl^kp when ^20 < 0, which provides the PI condition: 

^ <-^^(3cos^a- 1 ). (7) 

Eor g < 0 PI always leads to local collapses, resembling the 
situation in BECs with short-range attractive interactions [32, 
33] and in 3D homogeneous dipolar BECs [34, 35]. We hence 
focus below on BEC with repulsive contact interactions, g > 
0, where PI arises due to the attractive part of the anisotropic 
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DDL As shown in Ref. [28], for a - Q and a - njl the post 
PI dynamics is characterized by the formation of a transient 
gas of bright solitons. The soli tons then eventually undergo 
inelastic collisions, fusing together to form larger ones. The 
larger solitons may survive against collapse if the EEC density 
is low enough to hold the 2D criteria, \ij. 2 d\ ^ 


vestigated in Sec. V. 

IV. TWO DIMENSIONAL BRIGHT SOLITONS WITH 
TILTED DIPOLES 

In this section, we analyze the stability of 2D solitons in a 
dipolar EEC with tilted dipoles, using a 3D variational solu¬ 
tion and numerical NLGPE solutions. 


The PI regions as a function of p and a with g > 0 (Pig. 
2a) provide the first estimation for the stability regions of 2D 
bright solitons. In Pig. 2(a), there are two different regions in 
the a-p plane satisfying the PI criteria, separated at the magic 
angle = 54.7 degrees (or 0.95 radians). Por a < the 
PI requires gd < 0 which demands inverting the anisotropic 
character of DDI via rotating fields. This makes the direct 
cross-over between the two PI regions impossible by contin¬ 
uously varying a. Note that, in both regions, a sufficiently 
large [6| may lead to local collapses. This introduces an upper 
cut-off for [0| for the creation of stable solitons, which will be 
estimated in Sec. IV using variational calculations. Note from 
Fig. 2(a) that a dipolar EEC may be driven into PI just by tilt¬ 
ing the dipoles from an initially stable conhguration, without 
changing the interaction parameters or trap geometry. As it 
can be easily accomplished by changing the orientation of the 
externally applied field, we propose this as an alternative sim¬ 
ple method to generate 2D bright solitons in dipolar EECs, 
combined with an adiabatic approach. This is numerically in¬ 


A. Gaussian ansatz 


We consider the following Gaussian ansatz; 
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The dimensionless variational parameters: L^, Ly and Ly pro¬ 
vide the Gaussian widths along x, y and z. Introducing this 
ansatz in the energy functional 
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where g - gl and gd - gdl The 

minimum of E(Ljc,Ly,Lz) provides the equilibrium widths 
{w°, of the soliton. The absence of this minimum re¬ 

sults in two distinct types of instability. If the repulsive part of 
the interactions dominates, the soliton expands without limits 
on the xy plane (w^y oo), whereas dominating attractive 
interactions lead to collapse (w^y^ —> 0). 


B. Stability analysis and properties of 2D solitons 

For a = 0, we recover the isotropic scenario of Ref. [4], 
where stable isotropic solitons demand 
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which requires gd < 0, i.e. the inversion of the DDI. For large 
values of g, the above criteria reduces to p < -3/87r. Any a > 
0 breaks polar symmetry, and hence 3^ vv°. We quantify the 
2D soliton anisotropy with the aspect ratio y = where 

{/, j] e { 2 c,y) with w® > such that y < 1. 

For P <0, increasing a from zero towards leads to soli¬ 
ton elongation along y, Wy > w°, since the DDI becomes more 
attractive along y than along x, as schematically shown in Fig. 
3a. In addition, both and w® increase monotonously with 
a since the overall attractive interaction is reduced. As shown 
in Fig. 3c, y shows a non-monotonous character with a mini¬ 
mal value. Near the expansion instability, which happens at a 
g- and grf-dependent angle ae (< q;„, and shown by hlled cir¬ 
cles in Fig. 3c-f), the soliton anisotropy diminishes. The latter 
occurs because close to both are very large, the inter¬ 
action energy becomes hence very small, and as a result the 
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anisotropy diminishes. For/i > 0, for a = 7r/2 the solitons are 
maximally anisotropic with > w^, for any value of g and gd 
(see Fig. 3b). As expected, when a decreases from njl, the 
anisotropy decreases monotonously, until instability against 
expansion occurs at a critical angle. Figures 3c (d) and 3e (f) 
show the results for different g (gd) values with a fixed gd (g). 
In both figures we show the results for y obtained from the 
variation ansatz discussed above and from the numerical sim¬ 
ulation of the 2D NLGPE, which are excellent agreement in 
basically for all cases. As expected, the anisotropy increases 
with growing \gdlgV The fact that the values of g and gd for 
with the 2D soliton is stable depend on the tilting angle a, al¬ 
lows for the observation of solitons in parameter regimes in 
which solitons are unstable either for a = 0 [4] or a - njl 
[5], hence easing the experimental realization of 2D solitons. 
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Figure 3. The figures (a) and (b) show the equilibrium configurations 
of the soliton for {'i) p < Q {Q < a < a„) and (ii) p > 0 {a„ < a < 
njl). The thick (red) arrow shows the orientation of the dipoles in 
the BEC. For the case (i) the soliton is more elongated along the y 
axis and hence the aspect ratio is taken as y = wj./w®. In figures (c) 
and (d) y as a function of a is shown for the case (i). The value of 
a at which each curves terminates (a = a^), shown by filled circles, 
gives the critical angle for the expansion instability. For case (ii) 
the soliton is more elongated along the x axis, and hence we define 
y = w®/vv®, the corresponding plots as a function of a are shown in 
(e) and (f). In this case for a < the solitons are unstable against 
expansion. 


C. Soliton gas formation after phonon instability 

As mentioned above, PI in 2D dipolar BECs may be fol¬ 
lowed by the formation of a transient gas of bright solitons, 
instead of the collapse characteristic of short-range interact¬ 
ing BECs [28]. The post-instability evolution and the dynam¬ 
ics of the emergent soliton gas depend crucially on the tilting 



Figure 4. The snapshots of the post-PI dynamics in a 2D dipolar 
homogeneous condensate. The figures (a) and (b) are for g = 12, 
P = 0.28 and a = 1.35 radians at times t = 13.5/tUj and t = 21 re¬ 
spectively. Similarly, the figures (c) and (d) are forg = 40,/I = -0.29 
and a = 0.6 radians at times t = 31/tu, and t = 90/tj, respectively. 
Figures (a) and (c) are the transient stripes patterns with the disloca¬ 
tion defects shown in ellipses. Figures (b) and (d) are the temporarily 
arranged, unstable, ordered state of the soliton gas, which eventually 
merge each other during collisions. 


angle a. Eigs. 4(a) and 4(b) show the post-instability density 
patterns for g - 12, p - 0.28 and a - 1.35 radians at times 
t - 13.5/tUj and t - 21 respectively. Note that the the first 
stages of the post-instability dynamics are characterized by 
the formation of a transient stripe pattern that present disloca¬ 
tion defects [36]. At these dislocations two stripes merge into 
one. Note that for the case considered the DDI is more attrac¬ 
tive along the x-axis, which results in stripes almost parallel 
to the X direction. The initial formation of stripes is character¬ 
istic of in-plane anisotropy resulting from a > 0. Stripes may 
be observed as well for < 0, as illustrated in Eigs. 4(c) and 
(d) for p - -0.29 and a - 0.6 radians at times? = 31 jco^ and 
t - 90/cUj, respectively. In this case, the DDI is more attrac¬ 
tive along the y direction, resulting in stripes along the same 
direction. As shown in Eig. 4(b) and (d), the density stripes 
eventually break down into anisotropic solitons with major 
axis along the direction in which the DDI is more attractive, 
resulting in an unstable ordered state of solitons. The solitons 
eventually attract each other, fuse together and may remain 
stable depending on its condensate density. We have observed 
that the dynamics of the soliton gas inherently slowed down 
for a ^ 0 compared to the isotropic case (a - 0). This is be¬ 
cause, in the latter case the attractive forces between the soli¬ 
tons are rotationally invariant, whereas in the former case it is 
dominant along one particular direction due to the anisotropy 
of DDI, which restricts the motion in other directions. 
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Figure 5. The eigen frequencies (a) and the square of the compo¬ 
nents of the eigen-vectors (b) of two low-lying xy modes for p = 0.36 
as a function of a. For a = njl, the soliton is maximally anisotropic 
and the modes exhibit a pure x (solid lines) and y (dashed lines) char¬ 
acter. As a reduces, the solitons become less anisotropic and near the 
expansion instability the modes turn into breathing and quadrupole 
like ones, akin to that of the isotropic case. 
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D. Variational Calculations : Low-lying excitations 


At this point we examine the lowest-lying modes of 2D 
dipolar solitons using a variational method [37, 38], where 
we use a time-dependent Gaussian as the trial wave function; 

W-70('))^ 

iff(x,y,z,t)^A{t) e (H) 

ri=x,y,z 


where tjq, w,j, a,j and are the time dependent varia¬ 
tional parameters and the normalization constant A{t) - 
I yjWxWyW^. We consider the soliton is static in the xy 
plane and assume xo(f) = 0 andyo(f) = 0. The above ansatz is 
then introduced in the Lagrangian density of a dipolar BEC: 
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The Lagrangian is then obtained by integrating over the whole 
space i.e., L - J (PrX- We then obtain the corresponding 
Euler-Lagrange equations of motion for the time-dependent 
variational parameters, see App. A. Note that the centre of 
mass motion (along the z axis) is decoupled from the inter¬ 
nal dynamics of the soliton. Below we analyze in detail the 
two in-plane {xy-) modes of the condensate, which lie low¬ 
est in the excitation spectrum. Eor a = 0 (isotropic case), 
those are the breathing and quadrupole modes [4]. When 
a + 0, the character of the in-plane modes change, and in 
particular as the anisotropy of the soliton increases the modes 
in the xy plane decouple into pure x and y modes. In order 
to visualize it, we plot the square of the components of the 
eigen-vectors b = bxX + byy of the two planar modes for a 
particular case with p - 0.36, see Eig. 5(b) and the corre¬ 
sponding eigen-frequencies (w) are shown in Eig. 5(a). Eor 
P - 0.36 and a - njl, the condensate is more elongated along 
the x-direction and hence the lowest mode (solid lines in Eig. 
5) accounts mostly the oscillation of condensate width along 
the x-axis. When a reduces or dipoles tilting out of the xy 
plane, the anisotropy of the soliton decreases and the modes 
become more breathing-like (the lowest one) and quadru pole¬ 
like modes of the isotropic case, Eig. 5(b). Near the expan¬ 
sion instability they almost turn into breathing and quadrupole 
modes. The eigen-frequencies of the low lying modes are also 
crucial for the adiabatic control or preparation of the dipolar 
soliton, as discussed in Sec. V, where we have taken a = a(t). 
In particular, the frequency of the lowest mode determines the 
rate at which a should be varied in time in order for the system 
to remain in the instantaneous ground state of the system. 


E. Experimentally Relevant atomic systems 

Now we discuss the stability regions of the soliton in the 
state-of-the-art experimentally realized dipolar condensates of 
Cr, Dy and Er. In particular, we focus on the anisotropic soli¬ 
tons with y6 > 0, as it is experimentally less challenging. We 
take (jj^ - 2n X IkHz and N - 10000, resulting in an effec¬ 
tive dipolar interaction strengths of gd - 22 (Cr), 170 (Er), 
334 (Dy) for different atoms. The resulting phase diagrams 
obtained from Gaussian-energy calculations are shown in Eig. 
6 as a function of the v-wave scattering length a and the tilting 
angle a. These phase diagrams are of considerable interest 
if one would like to tune the system from soliton unstable to 
soliton stable region either by tilting the dipoles orientation 
or tuning a using the Eeshbach resonance. Below we demon¬ 
strate the adiabatic preparation of the soliton in a Cr BEC, us¬ 
ing experimentally realistic parameters, by tuning the tilting 
angle in time in which the dipolar BEC undergoes the transi¬ 
tion from externally confined BEC to a self trapped one. 


V. ADIABATIC PREPARATION OF AN ANISOTROPIC 
BRIGHT SOLITON 

The dependence of the BEC stability and of the soliton 
properties on the tilting angle a opens interesting perspectives 
for the preparation of 2D solitons by modifying in real time 
the tilting angle. In contrast to the previous sections we con¬ 
sider at this point the presence of an external shallow trap on 
the xy plane with frequencies a>x,y = 27r x 10 Hz. We consider 
a Cr-BEC of A = 6000 atoms with = 2.71 x 1 kHz, initially 
prepared with a tilting angle o', = 1.22 radians and a = 9ao, 
such that under these conditions there is no self-trapping (i.e. 
no 2D soliton), and the only conhnement is provided by the 
xy trap. We may then keep the interaction parameters intact. 
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Figure 6. The stability regions (shaded) for the anisotropic solitons as a function of tilting angle a and the i-wave scattering length a (in units 
of Bohr radius ao) for Dysprosium, Erbium and Chromium condensates, with w. = In x IkHz and N = 10000. The dipolar interaction is 
calculated with the intrinsic magnetic moments of each of the atoms. For large values of a the soliton become unstable against expansion and 
for lower values it undergoes the collapse instability. 


and slowly tune a to a final value (a/ = 1.32 radians) within 
the soliton-stability region. The trap on the xy plane may be 
then removed. The latter should be done slow-enough such 
that one avoids the signihcant creation of excitations in the 
condensate that may affect the self-trapping. The numerical 
results using the real-time evolution of the 2D NLGPE are 
shown in Fig. 7, where we monitor the widths of the con¬ 
densate in the xy plane as function of time. The self-trapping 
of the condensate in the xy plane is evident from the periodic 
oscillation of the condensate widths, after the removal of the 
harmonic conhnement. Note that the created soliton presents 
a slight breathing motion on the xy plane and the numerical 
results using the 2DGPE for the condensate density at various 
instants are shown in Fig. 7 a-c. 

We can hnd different regimes of adiabaticity for the cre¬ 
ation of bright soliton in dipolar EEC based on the different 
time scales in the system [39]. Here, we employ the condition 
Tij T <Si Tfj, where - max{l/wB(a(f))), the equiva¬ 
lent of quantum mechanical linear adiabatic time scale, deter¬ 
mined by the inverse of the lowest collective (breathing) mode 
(i)B and the nonlinear time scale = max{l/|yU 2 D(Q'(f))l) de¬ 
pends on the instantaneous chemical potential Note 

that in our case we vary the tilting angle from a, to a/ in 
which the chemical potential changes from a positive to neg¬ 
ative value continuously through fj. 2 D - 0 (Tf^ - oo). Thus the 
adiabatic criteria simply becomes « T, which guarantee 
us that we do not signihcantly populate the excitations dur¬ 
ing the period T [40]. In the particular case considered here, 
~ 0.05s. The soliton preparation involves two steps, hrst 
in which we vary a and the second involves the removal of xy 
confinement, both done linearly in time. We varied a in 0.3s 
from a, to a/ and the trap is then removed in a duration of 
0.5s. 


VI. CONCLUSIONS 

In conclusion, we studied the physics of 2D bright solitons 
in dipolar condensates as a function of the orientation of the 
dipoles with respect to the soliton plane. As we showed, the 
tilting angle may enhance the experimental possibilities to ob¬ 
serve the self-trapped matter waves as well as provide a probe 




Figure 7. The numerical results of 2D GPE for the preparation of 
an anisotropic soliton starting from a weakly confined [ojx.y = 2n x 
lOHz) Cr BEC with N = 6000, a = 9ao and a = 1.22 radians, (a) The 
initial condensate density, (b) after Is and (c) after 3.5s. The tilting 
angle a(l) is varied linearly from a,- = 1.22 to Uf = 1.32 radians in 
0.3s. Then, the trap in the xy plane is removed in 0.5s. (d) The widths 
Wx and Wy of the soliton along x and y directions respectively as a 
function of time. The periodic oscillation of the widths for t > 0.8s 
indicates the self-trapping of the condensate. 


to tune the anisotropy of the soliton. In addition, it can drive 
the dipolar BEC into PI without altering its interaction param¬ 
eters or trap geometry. In 2D, the post PI dynamics is always 
characterized by a transient stripe formation and eventual for¬ 
mation of an unstable soliton gas. Finally we have demon¬ 
strated how to prepare the 2D soliton adiabatically by tuning 
the tilting angle in the case of a Cr BEC. 
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The above equations describe the motion of a particle with 
coordinates in an effective potential 


Appendix A: Lagrangian for a dipolar condensate 


The Lagrangian density for a dipolar condensate is given in 
Eq. 12, and using the Gaussian time-dependent trial function 
we obtain the Lagrangian L - J cPr£.: 




ri=x,y,z 

, g 


.Priwl 

2 2m 


v2w2 


■wy. 


1 2^1 


1 


^J2^ 4;rw^Wy Wj 


+ E(w„), 


(Al) 


where 


f n e-^ (A2) 

^(2n) J ^ t. 


with Vd(k) the Eourier transform the dipole-dipole potential. 
Then, the equations of motion for the condensate widths are 
obtained as 


mWx = -r -H 

mWy 


fp- 

iw 


mwy - -T -H 

mWy 


''z - 3 

mw. 


8 


(A3) 

{27T)'^/^wlWyWz 

dwx 

8 


(A4) 

(27Ty/^WxWyWz 

1 

^1 

8 

2 2 

-2^. (A5) 
dwz 

{27T)'^^^WxWyWl 



= I; E 4 * + 

ri V 


Cln'fl'^WxWyWz 


(A6) 


Once the equilibrium widths of the condensate are obtained 
by minimizing the effective potential (or equivalently from the 
Gaussian energy calculations in Section. IV A), the low lying 
excitations are obtained by diagonalizing the Hessian matrix 
of U. Also, note that the centre of mass motion of the soliton 
along the z axis is de-coupled from the internal dynamics and 
is governed by the equation 


zo 


-OJ^ZO- 


(A7) 
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